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Context. The presence of electric currents in the atmospheres of magnetic chemically peculiar (mCP) stars could bring important 
' theoretical constrains about the nature and evolution of magnetic field in these stars. The Lorentz force, which results from the inter- 

action between the magnetic field and the induced currents, modifies the atmospheric structure and induces characteristic rotational 
variability of pressure-sensitive spectroscopic features, that can be analysed using phase-resolved spectroscopic observations. 
Aims. In this work we continue the presentation of results of the magnetic pressure studies in mCP stars focusing on the high- 
resolution spectroscopic observations of Bp star 56 Ari. We have detected a significant variability of the Ha, H/3, and Hy spectral 
lines during full rotation cycle of the star. Then these observations are interpreted in the framework of the model atmosphere analysis, 
which accounts for the Lorentz force effects. 

Methods. We used the LLmodels stellar model atmosphere code for the calculation of the magnetic pressure effects in the atmosphere 
of 56 Ari taking into account realistic chemistry of the star and accurate computations of the microscopic plasma properties. The 
SYNTH3 code was employed to simulate phase-resolved variability of Balmer lines. 

Results. We demonstrate that the model with the outward-directed Lorentz force in the dipole+quadrupole configuration is likely to 
J^l . reproduce the observed hydrogen lines variation. These results present strong evidences for the presence of non-zero global electric 

■ currents in the atmosphere of this early-type magnetic star. 

■ Key words, stars: chemically peculiar - stars: magnetic fields - stars: atmospheres - stars: individual: 56 Ari 
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^ — 1 ' 1. Introduction ILehmann et al.ll2007l) . At the same time, some of the magnetic 

stars demonstrate the characteristic shape of the Balmer line 

O ■ The atmospheres of magnetic chemically peculiar (mCP) stars variability that can not be simply descri bed by temperature or 

L; ; display the presence of global magnetic fields ranging in strength abundance effects. For example, Krolj (TT^89h showed that at 

. , _h , from a few hundred G up to several tens of kG (| Landstreel]2001|) , leafjt part of the variability detected in several mCP stars can 

^ ■ with the global configurations well represented by dipolar or be attributed to the pressure effects, indicating the presence of a 

H I low-order multipolar components (|Bagnulo et al.||2002D and that non . zero Lorentz force in their atmospheres. 

5t i are likely stable during significant time intervals. The stability 

of the atmospheres against strong convective motions and the Different atmosphere models with the Lore nz force were 
existence of large scale magnetic fields provide a unique condi- consi dered by several authors (see review in [Valyavin et ail 
tions for the study of secular evolution of global cosmic mag- 2004 ). In th i s stu dy we follow approaches given by 
netic fields and other dynamical processes which may take place IValvavin et ail d2004l) considering the problem in terms of in- 
in the magnetized plasma. In particular, the slow variation of the duced atmospheric electric currents interacting with magnetic 
field geometry and strength changes the pressure-force balance fields. The authors predicted that the amplitude of the variations 
in the atmosphere via the induced Lorentz force, that makes it in hydrogen Balmer lines seen in real stars can be described 
possible to detect it observationally and establish a number of if ne assumes strong electric currents fl owing in upper atmo- 
important constraints on the plausible scenarios of the magnetic spheric layers of these objects. Later on, IShulvak et all (120071) 
field evolution in early-type stars. (Paper I hereafter) extended and improved this model to more 
Among the known characteristics of mCP stars, the vari- complicated field geometries and accurate treatment of mag- 
ability of hydrogen Balmer lines is poorly understood. For netized plasma properties. Their first direct implementation of 
some of stars it can be connected with inhomogeneous surface new model atmospheres to the analysis of the hydrogen spectra 
distribution of chemical elements, possible temperature vari- of one of the brightest mCP star 8 Aur showed that their rota- 
ations and/or stellar rotation (see, for example, discussion in tional modulation can be induced by the Lorentz force effect, 

which is not directly connected to the temperature or abundance 

Send offprint requests to: D. Shulyak, variation across the stellar surface. The precise analysis of the 

e-mail: denis.shulyak@gmail.com longitudinal magnetic field variation and subsequent modeling 
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of the magnetic pressure for every observed rotational phase of 
the star allowed us to constrain the magnitude and the direction 
of the Lorentz force. In particular, the outward-directed Lorentz 
force (i.e. directed outside the stellar interior along radius) was 
found to provide the best fit to observations, in combination with 
rather strong induced e.m.f. (electro-magnetic force) of about 
1 x 10" CGS units which was found to play an important role 
in the overall hydrostatic structure of the stellar atmosphere. 

The knowledge of the direction and the strength of the 
Lorentz force is important for the understanding of physical 
mechanisms that are responsible for such strong surface currents 
and their interaction with global, large-scale magnetic fields (see 
discussion in Paper I for more details). Taking this in to account 
and following the pioneering work bv lKroif (fT989h . we have 
initiated a new spectroscopic search of hydroge n line variabil- 
ity in a number of magnetic main-sequence stars dValvavin et al.l 
120051) . 

In this paper we present the phase-resolved high-resolution 
observations of one of the weak-field (|(B Z )| < 500 G) mCP star 
56 Ari (HD 19832). We detected significant variation of the 
Balmer line profiles and interpreted it in terms of the non-force- 
free magnetic field configuration. 

The overview of observations will be presented in the next 
section. Then, in Sect. 3 we will give a short description of 
the model used to simulate the effects of the magnetic pressure. 
Main results will be summarized in Sect. 4 with conclusions and 
discussion given in Sect. 5 and 6 respectively. 



2. Observations 

Observations of 56 Ari were carried out with the BOES echelle 
spectrograph installed at the 1.8 m telescope of the Korean 
Astronomy and Space Science Institute. The spectrograp h and 
ob servational proced ures are described bv lKim et al. (2000) and 
bv lKim et al.l (|2007). The instrument is a moderate-beam, fibre- 
fed high-resolution spectrometer which incorporates 3 STU 
Polymicro fibres of 300, 200, and 80 ymx core diameter (cor- 
responding spectral resolutions are A/AA = 30000, 45 000, 
and 90 000 respectively). The medium resolution mode was em- 
ployed in the present study. Working wavelength range is from 
3500 A to 10000A. 

Seventeen spectra of the star were recorded in the course of 
10 observing nights from 2004 to 2006. Typical exposure times 
of a few minutes allowed to achieve S /N ~ 250-300. Table Q] 
gives an overview of our observat ions. Throughout this study 
we implement ephemeris derived by Adelm an et alj d2001l) using 
linear changing period model: 



Table 1. Observations of 56 Ari. 



JD = 2434322.354 + 



0.7278883 



1 +0.7278883-5 -(f-fo)' 



(1) 



where t = 2434322.354 and S = -1.35 • 10~ 9 . 

Details related to spectral data reduction and processing, 
also study of the spectrograph's stability are presented in 
IShulvaketal.1 (120071) . and we don't describe it here. Accuracy of 
the continuum normalization around Balmer lines is estimated to 
be approximately 0.2 - 0.3%. 

3. Model 

3.1. General equations and approximations 

In this section we follow the approach and methods outlined pre- 
viously in Paper I. However, for the sake of explanation, we find 



No. 


JD 


Rotation Phase 


1 


2453250.2191 


0.321 


2 


2453250.3017 


0.435 


3 


2453251.2284 


0.708 


4 


2453251.3359 


0.855 


5 


2453306.1282 


0.128 


6 


2453306.2942 


0.356 


7 


2453308.0681 


0.793 


8 


2453309.0295 


0.114 


9 


2453309.1620 


0.296 


1U 


T/ici^t noon 
245.3 O66.0S5U 




11 


2453671.0690 


0.478 


12 


2453671.0860 


0.502 


13 


2453759.0098 


0.290 


14 


2453760.0284 


0.689 


15 


2453760.9714 


0.985 


16 


2453762.0026 


0.402 


17 


2453762.9682 


0.728 



it useful to state here some of the general assumptions used in 
the modeling procedure: 

1. The stellar surface magnetic field is axisymmetric and is 
dominated by dipolar or dipole+quadrupolar component in 
all atmospheric layers. 

2. The induced e.m.f. has only an azimuthal component, sim- 
ilar to that described by Wrubel ( 1952), who considered de- 
cay of the global stellar magnetic field. In this case the dis- 
tribution of the surface electric currents can be expressed by 
the Legendre polynomials P l n {p), where n — 1 for dipole, 
n — 2 for quadrupole, etc., and [i = cos is the cosine of the 
co-latitude angle 9 which is counted in the coordinate system 
connected to the symmetry axis of the magnetic field. 

3. The atmospheric layers are assumed to be in static equilib- 
rium and no horizontal motions are present. 

4. Stellar rotation, Hall's currents, ambipolar diffusion and 
other dynamical processes are neglected. 

Taking these approximations into account and using 
Maxwell equation for field vectors and Ohm's law, one can write 
the hydrostatic equation in the form 



dPn 



dr 



= -pg ± -A ± J] CnPfo) J] ^ = ~ 



•Pgeff- 



(2) 



Obtaining this equation we used the superposition principle for 
field vectors and the solution of Maxw ell equa t ions f or each of 
the multipolar components following Wrubel (1952). We also 
suppose that E _L B. Here c„ represents the effective electric 
field generated by the n-th magnetic field component at the stel- 
lar magnetic equator and B B is the horizontal field component. 
The signs "+" and "-" refer to the outward- and inward-directed 
Lorentz forces, respectively. 

We note that the values of c„ are free parameters to be found 
by using our model. These values represent the fundamental 
characteristics that can be used for building self-consistent mod- 
els of the global stellar magnetic field geometry and its evolu- 
tion. Thus, an indirect measurement of these parameters via the 
study of the Lorentz force is of fundamental importance for un- 
derstanding the stellar magnetism. 

Calculation of the electric conductivity A ± is carried out us- 
ing the Lorentz collision model where only binary collisions be- 
tween particles are allowed which is a good approximation for a 
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Table 2. Abundances (in log N^/Num) of 56 Ari, used for deter- 
mination of model atmosphere parameters. 





He 


Mg 


Al 


Si 


Fe 


56 Ari 
Sun 


-2.10 
-1.10 


-5.51 
-4.51 


-6.17 
-5.67 


-3.53 
-4.53 


-4.09 
-4.59 



low density stellar atmosphere plasma. The detailed description 
and basic relationships of this approach are given in Paper I. 

From the Eq. (O it is seen that, in the presence of electric cur- 
rents and magnetic field, the rotation of a star can produce phase- 
dependent Lorentz-force term (due to a variation of B^ and A±), 
which will in turn modify the hydrostatic structure of the at- 
mosphere, manifesting itself as a variation of pressure-sensitive 
spectral lines. 

3.2. Model atmospheres with Lorentz force 

Our calculations were carried out wit h the stellar model a tmo- 
sphere code LLmodels developed by Shulvaketal. (2004). At 
each iteration the code calculates electric conductivity in all at- 
mospheric layers using all available charged and neutral plasma 
particles. The conductivity is then used to evaluate the magnetic 
contribution to the magnetic gravity and to execute temperature 
and mass correction procedure. 

The input parameters for the calculation of magnetic pres- 
sure are the direction of the Lorentz force (inward- or outward- 
directed), e.m.f. at the stellar equator, mean surface magnetic 
field modulus (B), and the product of the two sums in Eq. (O 
containing contribution from all considered multipolar compo- 
nents for every single rotational phase of the star. 

As can be seen from Eq. (0, there is some critical value of 
c„ that may produce unstable solution in the case of the outward- 
directed Lorentz force. Such models cannot be considered in 
the hydrostatic equilibrium approximation introduced above and 
were deemed non-physical in our calculations. Thus, for each set 
of models, we restricted c„ values to ensure static equilibrium. 

The atomic line lis t was extracted from the VALD database 
dPiskunov et all 119951; iKupka et all 1 19991) . including all lines 
originating from the predicted and observed energy levels. This 
line list was used as input for the lines opacity calculation in the 
LLmodels code. 

4. Numerical Results 

4.1. Model atmosphere parameters of 56 Ari 

The model atmosphere parameters, log(g) and T e s, were deter- 
mined using theoretical fit of the hydrogen B aimer lines and 
spectral energy distribution. It was constructed by combining 
the average of the optical spectrophotometric scans obtained by 
I Adelmarjdl 983b and low dispersion UV spectrograms from the 
IUE INES3database. For modeling the hydrogen HB and Hy line 
profiles we used the mean spectrum of 56 Ari averaged other all 
17 observed rotational phases. The projected rota t ional velocity 
y sin; = 160 km s _1 was taken from lHatzesI (1 1 9931) . iRy a bchikova 
(2003) used the same value in a more recent Doppler imag- 
ing study. Individual abundances of several chemical elements, 
listed in Table|2] were determined as described below (Sect. l4~3l l. 

Synthetic Balm er line profiles we re calculated using the 
SYNTH3 program (lKochukhovl l2007). which incorporates re- 

1 |http : // ines . ts . astro . itj] 



cent improvements in the treatment of the hydrogen line opac- 
ity (jBarklem et al. 2000). The stellar energy distribution and 
Balmer lines are approximated best with the following param- 
eters: r eff = 12 800 + 300 K, log(g) = 4.0 + 0.05. Note that 
such high accuracy of the determined parameters is just an in- 
ternal accuracy obtained from our technique, which we used to 
fit the data. Real parameters may be slightly different from the 
obtained ones due to various systematic error sources, but this 
does not play a significant role in our study. 

Comparisons of the observations and model predictions 
are presented in Figs. Q] and [2] We transformed Adelman's 
spectrophotometr i c obse rvations to the absolute units following 
Lipski & Stepieri (2008). Since the absolute calibration of IUE 
fluxes around their red end have substantial uncertainties (see 
iLipski & Stepienll2008l) . we scaled the IUE fluxes by * 10% to 
match the Adelman's data in the near-UV region. This correction 
is comparable to the offset b etween alternative flux calibrations 
suggested for the IUE data dGarcia-Gil et al.l [20051) . Applying 
the correction ensured that the IUE spectra are smooth contin- 
uation of the optical spectrophotometry. Then the theoretical 
fluxes can be adjusted to fit the combined set of observations. 
Lipski & Stepieh (2008) also noted the discrepancy between ob- 
servations in visual and UV, but did not do any attempts to cor- 
rect it. This is why their final effective temperature of 56 Ari was 
found to be T e ff = 12 250 K resulting from the fact that ignor- 
ing the obviously spurious offset between the observed datasets 
one needs a lower r e (j to fit the IUE fluxes simultaneously with 
Adelman's fluxes redward of the Balmer jump. As an exam- 
ple we also show in Fig. Q] theoretical flux obtained from the 
T e ff = 12 300 K, log(g) = 3.9 model which fits reasonably well 
the hydrogen Balmer lines (see Fig. [2j, but fails to reproduce 
either Adelman's or IUE data. 

Photometric observations in the Stromgren and UBV sys- 
tems also point to the higher T e ff and log(g) of 56 Ari. For 
instance, comp aring theoretically computed col or-indices with 
observations of lHauck & Mermilliodl (U998) and iNicoletid 19781) 
we find that observed photometric parameters (b - y - -0.052, 
ci = 0.55, B — V = -0.12, U - B = -0.42) are best fitted 
with r eff = 12 800 K, log(g) = 4.0 model (b - y = -0.052, 
ci = 0.60, B - V = -0.12, U - B = -0.43) rather than more 
cooler r eff = 12 300 K, log(g) = 3.9 model (b - y = -0.048, 
ci = 0.67, B -V - -0.12, U - B = -0.38). 

Finally, recent stu dies by Kochuk hov et al.l d2005l) and 
iKhan & S hulyak (2006) showed that the effects of Zeeman split- 
ting and polarized radiative transfer on the model atmosphere 
structure and shapes of hydrogen line profiles are less than 0.1% 
for magnetic field intensities around 1 kG, so they can be safely 
neglected in the present investigation. 

4.2. Magnetic-field geometry 

To calculate the Lorentz force effects, it is essential to specify 
the magnetic-field geometry (see Eq. (O). For this purpose we 
made use of the longitudinal magne tic field measurements ob- 
tained by iBorra & Landstreel d!980l) using HB photopolarimet- 
ric technique. The authors observed a smooth single-wave (B z ) 
variation with rotation phase and concluded that it is probably 
caused by a dipole inclined to the rotation axis of the star. 

Generalizing this work, we have approximated the magnetic- 
field topology of 56 Ari by a combination of the dipole and ax- 
isymmetric quadrupole magnetic components. We also assumed 
that the symmetry axes of the dipole and quadrupole magnetic 
fields are parallel. Thus, the magnetic model parameters include 
the polar strength of the dipolar component B^, relative con- 
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1000 2000 3000 4000 5000 6000 7000 

A, [A] 

Fig. 1. Comparison of the observed and computed spectral energy distributions of 56 Ari. Theoretical models correspond to T e g = 
12 300 K, log(g) = 3.9 and T eS = 12 800 K, log(g) = 4.0. The model fluxes have been convolved with an FWHM = 10 A Gaussian 
kernel for a better view. 




AX, [A] 

Fig. 3. Comparison of the longitudinal field observations of 
Fig. 2. Comparison of the observed and computed Ha, Hp, and 56 ^ (symbo i s) and the mo del curves for i = 60° and two 
Hy line profiles. extreme values of the B q /Bd parameter (lines). 



tribution of the quadrupole field fi q /fid, magnetic obliquity B, 
and inclination angle i of the stellar rotation axis with respect to 
the line of sight. The last parameter can be estimated from the 
usual oblique rotator relation connecting stellar radius, rotation 
period, and vsin/. Employing the recentl y revised Hipparco s 
parallax of 56 Ari, n - 6.49 + 0.76 mas (Ivan Leeuwenll2007l) . 
r e ff = 128 00 K and bolome t ric cor rection BC = -0.74 deter- 
mined by iLipski & Stepienl d2008l) . we found a stellar radius 
R = 2.8 + 0.4 R e and inclination angle i = 57 + 13 ° . How ever, 
the previous Dopp ler imaging studies of iHatzesi d 19931) and 
Rvabchikova ( 2003) favoured the value of i = 70°. Consequently 
we decided to explore the model parameters for the entire range 
of i = 50-70°. 



The remaining free parameters of the magnetic field model 
were dete rmined with the least-squa re fit of the observed {B z ) 
variation dBorra & Landstreeil Il980h . Compared to Aur in- 
vestigated in Paper I, the longitudinal magnetic field curve of 
56 Ari is poorly defined due to large observational errors and 
a relatively small number of measurements. In this situation 
the acceptable solutions for the magnetic field geometry fall in 
a rather wide range: from B q /Bd = 0.5 to B q /B,\ = -5 for 
i = 50°, to B q /B d = -4.5 for i = 60°, and to B q /B d = -4 
for i = 70° (Fig. |3). The corresponding dipolar field strength 
range is B& - 1.3-1.8 kG and magnetic obliquity is B — 70-90°. 
The lowest x\ = X 1 l v was found between fiq/Ba = [—1-5, -3] 
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with the average valued « 0.66 while xl ~ 0.71 for other ge- 
ometries: with such a small difference we conclude that all the 
considered magnetic field geometries are equally possible. 

4.3. Effects of horizontal abundance distribution 

We used spectrum synthesis calculations to access chemical 
properties of the atmosphere of 56 Ari. Abundances of He, Al, 
Mg, Si, and Fe were determined by fitting SYNTH3 spectra 
to the average observations of the star. Results of this analy- 
sis, summarized in Table [2] indicate that He, Al and Mg are 
deficient in the atmosphere of 56 Ari with respect to the solar 
chemical composition. Fe is moderately overabundant while Si 
is strongly enhanced. Abundance of Cr cannot be determined re- 
liably, but since no prominent Cr n lines are present in the spec- 
trum of 56 Ari, we concluded that its overabundance is most 
likely smaller than that of Fe. 

In addition to the mean abundances analysis we interpreted 
variation of the pha se-resolved spectra usin g the Doppler imag- 
ing (DI) technique (Kochukhov et al. 2004). Details of this work 
will be presented in a separate publication. In the present paper 
we are only concerned with the possible effect of the horizontal 
abundance variations on the hydrogen line profiles. Effects of the 
possible vertical stratitification of chemical elemets are ignored 
because for the purpose of our study they will not differ much 
from the effects of horizontal abundance inhomogeneities. By 
deriving surface-resolved abundances we effectively account for 
the line opacity variation which would be introduced by chemi- 
cal stratification. 

Among the elements mentioned above only Mg and Si show 
strong line profile variations. Therefore, we reconstructed abun- 
dance maps of Si using the red doublet Sin 6347, 6371 A and 
Mg using the 4481 A line. The modeling of the latter region 
also included the He 1 4471 A line, which allowed us to estimate 
potential influence of the He abundance variation on the model 
atmosphere structure and hydrogen line profile behaviour. 

Analysis of Si, He, and Mg demonstrated that horizontal 
abundance inhomogeneities give a negligible contribution to the 
hydrogen line profile variation. We have computed model at- 
mospheres with the surface-averaged abundances for each ro- 
tational phase of the star. Figure [4] shows that the resulting stan- 
dard deviation of the synthetic profiles is minute compared to the 
variation observed in the line wings of Ha, HB, and Hy. 

Much larger abundance gradients of the iron peak elements 
are required to induce a noticeable modulation of the hydro- 
gen line profiles during rotation cycle. The dashed line in Fig. [4] 
shows results of the numerical experiment where we have as- 
sumed that Fe has the same horizontal distribution as Si but with 
10 times higher contrast. This calculation contradicts the actual 
observations since the Fe lines in 56 Ari are varying weakly and 
differently from Si, but it represents a useful illustration of the 
impact of chemical inhomogeneities. One can see that the metal 
abundance spots lead to the variation of the depth of hydrogen 
line cores, while in 56 Ari we see chang es in the line wings. 
Thus, following Paper I and iKrollI £[989), we are led to con- 
clude that the chemical spots can not contribute to the observed 
Balmer line wing variations. 

4.4. The Lorentz force 

To predict the phase-resolved variability in hydrogen line pro- 
files both the inward- and outward-directed Lorentz forces 
are examined through the model atmosphere calculations. The 
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Fig. 4. The standard deviation cr at the Hot, HB, and Hy lines. 
Observations are shown by solid lines. Dotted line - variability 
due to inhomogeneous element distributions in 56 Ari. The ad- 
ditional effect of the hypothetical extreme horizontal Fe gradient 
is shown by the dashed line (see text for explanations). 



actual magnetic input parameters of computations with the 
LLmodels code include the sign of the Lorentz force, magnetic 
field modulus B, and the product of two sums 

2 c n P\ Z Bf (see 

Eq. [2). We take the last two parameters to be disk-averaged at 
the individual rotation phases incorporating them in ID stellar 
atmosphere code. The corresponding phase curves of the mag- 
netic parameters relevant for our calculations are illustrated in 
Fig.i 

Taking into account the large variety of possible solutions for 
the surface magnetic field geometry (see Sect. |4.2| > we calculated 
a number of model grids with the outward- and inward-directed 
Lorentz force and B q /Bd ranging from 0.5 to -4.5 with a step 
of 0.5. 

The calculations showed that, in order to reproduce the am- 
plitudes of the observed standard deviations due to the profile 
variations, the effective electric field should be in the range 
c\ = lxl0~ 10 -5 xl0~ 10 CGS in the case of the inward-directed 
Lorentz force and cj = 1.2x 10" 11 - 1 x 10" 10 CGS in the case of 
the outward-directed Lorentz force. Under the assumption of a 
purely dipolar configuration these values are: c\ = 5xl0~ 10 CGS 
for inward-directed and c\ = 7.5 x 10~ n CGS for outward- 
directed Lorentz forces. 

Not for all models considered in our study it was possi- 
ble to fit the amplitude of the observed standard deviation with 
the inward-directed Lorentz force. This is true for models with 
Bq/Ba > -1, where the changes in longitudinal magnetic field 
and magnetic field modulus result in very narrow phase-resolved 
variations in magnetic force term. As an example, the left panel 
of Fig. illustrates the run of the g e ff with the Rosseland op- 
tical depth in the atmosphere of 56 Ari (Z? q /Z?d = -2.5) for 
inward- and outward directed Lorentz forces computed under 
different assumptions about induced electric field. An increase 
of the e.m.f. value by a factor of two considerably changes the 
amplitude of the effective gravity, but the computed standard de- 
viation does not change much (see lower left panel of Fig. [5}. 
This is why even a large increase of e.m.f. does not reveal it- 
self in a standard deviation plot. Obviously, the situation may 
change once a more complex geometry of the magnetic field is 
introduced, however it is connected with the introduction of ad- 
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hward-directed Lorentz force, B q /B =-2.5, c ; /c, = 7 Inward-directed Lorentz force, B Q /B d =-0.5, c,/c, = 0.5 Inward-directed Lorentz force, B q /B d =0.0, c,/c, = 0.0 




-30 -10 10 20 -20 -10 10 20 -20 -10 10 20 

AX (A} AX (A) AX (A) 



Fig. 5. The effective acceleration as a function of the Rosseland optical depth for different rotation phases calculated for several 
magnetic field configurations and induced e.m.f. The resulting standard deviations around Hy line are shown in the bottom panel. 



ditional free parameters making the fitting procedure ambiguous. 
In contrast, the outward-directed magnetic force seems to have 
a larger impact on the model pressure structure: varying e.m.f. 
value from 1 x 1CT 11 CGS to 1.2 x 10" 11 CGS changes the am- 
plitude of standard deviation by about a factor of two or more. 

The amplitude of the standard deviation around hydrogen 
lines strongly depends on the magnitude of the inward-directed 
Lorentz force for different magnetic field geometries. For in- 
stance, the middle panel of Fig. [5] illustrates standard deviations 
for the model with B q /Bd = -0.5: changing induced e.m.f. by 
a factor of two considerably increases the amplitude of the stan- 
dard deviation. Similarly, the right panel of Fig.[5]shows the pre- 
dicted variations for the purely dipolar model. 

4.5. Comparison with the observations 

In the following we compare the residual theoretical and ob- 
served Balmer lines with model predictions based on different 
assumptions about the magnetic field geometry and the direction 
of the Lorentz force. The residuals are obtained by subtracting a 
spectrum at a reference phase (<p = 0.1 14 where the Balmer pro- 
files have the largest widths) from all the other spectra. Figure [7] 
illustrates residual Ha, H/3, and Hy line profiles for each of the 
observed rotation phases. The positive sign of the residuals im- 
plies that the lines at the current phase are narrower than those 
obtained at the reference phase. It is seen that the characteristic 
behaviour of hydrogen lines demonstrates a single wave varia- 
tion with the most noticeable effect at phases between cp - 0.402 
and <p - 0.855. The effect is also seen in the red wings of lines 
at (p — 0.985, however it is smeared out in the blue wing. We 
emphasize that this systematic asymmetry, with the blue wing 



lying below the red one, is observed for all three studied Balmer 
lines. Inaccuracy of the spectrum processing can be one of the 
reasons for this effect. Howe ver, our data reduct ion is identical 
to the analysis of 9 Aur (see Shu lyak et al.ll2007l) . which shows 
no such asymmetry. Thus, we suspect that the asymmetry may 
also have physical origin due to non-stationary, magnetically- 
channeled stellar wind from the surface of 56 Ari. Very fast ro- 
tation and a relatively high temperature of this star make it plau- 
sible that its wind produces the variable, obscured P-Cyg fea- 
ture distorting blue wings of the Balmer lines. The presence of a 
weak feature approximately 2. 5 A blueward the line center seen 
for all three hydrogen lines could also be an argument for the 
wind (for Hy line the presence of this feature could alternatively 
be explained by absorption in Ti n and Fe n lines, however in case 
of H/3 there are no spectral lines that could contribute to this fea- 
ture). Nevetheless, despite these problems introduced by an un- 
known physical process, the characteristic shape of the Lorentz 
force induced variability is clearly seen in the hydrogen lines of 
56 Ari, making it possible to perform the analysis in the frame- 
work of our modeling approach. 

None of the tested theoretical models appeared to fit the ob- 
served variability in all phases, but still generally describes the 
data more or less reasonably well for more than half of them 
(see, for example, Model 1 or Model 4). Varying the magnetic 
field geometry and the ratio of induced dipolar and quadrupolar 
equatorial e.m.f. 's (C2/C1), we could achieve a good agreement 
for a certain phase interval only: either it was possible to fit the 
observations around phase <p = 0.5 or for other phases only. As 
an example, in Fig. [7] we plot some of the theoretical predictions 
for models that provide a more or less reasonable fit to the obser- 
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vations. In particular, models with the outward-directed Lorentz 
force are shown for the following configurations: 

B q /B d = -2.5, ci = 1.2XKT 11 CGS, c 2 /c { = 7.0 (Model 1), 
B q /B d = -4.5, a = 1.2xl0 -11 CGS, cz/c\ = 7.0 (Model 2), 
B q /B d = +0.5, ci =9.0x10-" CGS, c 2 /c l = 0.1 (Model 3). 

Also, the model with the inward-directed Lorentz force is pre- 
sented: 

B q /B d = -1.0, ci = 5 x 10- 10 CGS, c 2 /c\ = 2 (Model 4) 

For all plotted models the inclination angle i = 50° was assumed. 
Taking i - 70° does not change much the disc-integrated param- 
eters of the magnetic field and thus leads to essentially the same 
picture of the hydrogen lines variation (see below). 

Models 1 and 2 have the same parameters except the strength 
of quadrupolar magnetic field component. They produce almost 
the same fit to the observed variations, and in the same man- 
ner are not able to fit phases at <p — 0.636 and above (Although 
Model 2 gives systematically a little bit better fit there). Note 
that models with B q /B d < —1.5 give the same kind of the 
fit, but we do not plot them here to avoid overcrowding the 
plot. At the same time, Model 3 seem to be a preferable one 
for these phases, however it fails to fit observations at phases 
<p - 0.435,0.478,0.502, and gives enormously high effect at 
phases <p = 0.321 and down to zero. This is true also for the 
Model 4 with inward-directed Lorentz force. This model fits rea- 
sonably well such phases as <p — 0.689, 0.708, but yields the 
line wings that are generally too wide comparing with observed 
ones. Thus, of the two possible directions of the Lorentz force 
in our model we consider an outward-directed Lorentz force as 
the more reasonable choice to describe observations of 56 Ari. 
Due to problems with telluric lines the continuum normalization 
around Ha line is substantially inaccurate comparing to other 
lines and it is not possible to distinguish between different mod- 
els there. 

Testing models with different magnetic field parameters we 
tried to find those that predict a single-wave variation of the mag- 
netic force term over the rotation cycle, as indicated by observa- 
tions. Moreover, its run is likely to have a wide plato around 
<p = 0.5 and drop rapidly close to <p = and <p = 1 to fit obser- 
vations (see Fig. |7]i. By varying the parameters B q /B d and C2/C1 
(with the fixed i and B) we succeeded to find sets of parameters 
that give this kind of plato, but in all cases it appears to be not 
as wide as it is needed to fit observations in all phases. This is 
illustrated in Fig.|6]where we plot magnetic parameters used for 
the Lorentz force calculation in some of the models mentioned 
above as a function of the rotation phase. The right panel in this 
figure illustrates predictions for the purely dipolar model. It is 
also seen that the inclination angle i does not play critical role in 
the present investigation: models with different i would give the 
same phase-resolved variation in hydrogen lines and any ampli- 
tude difference between them can be adjusted by a proper choice 
of a c\ parameter. 

In this investigation we focused analysis on the hydrogen 
lines. It appears that they are most sensitive to the pressure 
changes introduced by the Lorentz force. As for metal lines, 
none of the strong Si n lines visible in the spectrum of 56 Ari, 
A 5055.98 A, 5466.48 A, 5466.89 A, 6347.11 A or 6371.37 A, 
exhibits significant variation due to a non-zero Lorentz force. We 
have tested this by using the magnetic parameters of Model 1 
(which has the largest amplitude of the g e ff variation, see Fig. [5j 
and recomputing spectrum models for every phase with the mean 
abundances. We find no detectable changes in the line wings and 



less than 1% difference in the line cores between models with 
and without Lorentz force. This difference is likely to be due to 
the differences in the temperature distribution of these two mod- 
els. No visible phase-dependent changes can be seen in the spec- 
tra corresponding to the models with Lorentz force. These results 
lead us to the conclusion that, due to their high pressure sensi- 
tivity in the predominantly ionized plasma of the atmosphere of 
such relatively hot star, only the hydrogen lines are useful indi- 
cators for the magnetic pressure effects. 

Similarly, we find no evidence for the influence of phase- 
dependent pressure effects on the stellar spectral energy distribu- 
tion. The maximum difference between models with and without 
Lorentz force is less than 2% in Balmer continuum. This corre- 
sponds to « 0.01 mag difference in c\ color-index and even less 
for other Stromgren parameters. Thus, variation seen, for exam- 
ple, in the ph ase-resolv e d sp ectrophotometric scans of 56 Ari 
published by Adelman (1983) could not be attributed to the 
Lorentz force effects but are produced by inhomogeneous abun- 
dances and/or other mechanisms. 

We note that a strong decrease of g^ evident in Fig. (up 
to 1 dex around log(rR OSS ) = 0) compared to the non-magnetic 
case leads to a relatively small difference in the observed pa- 
rameters due to a) the fact that this decrease does not affect the 
entire stellar atmosphere and b) non-local nature of the hydro- 
static equation in the presence of depth-dependent g e $. The later 
implies that, for example, one order of magnitude increase of the 
magnetic gravity results only in three times lower gas pressure 
for the outward-directed Lorentz force, which is too small to sig- 
nificantly change the opacity coefficient and influence the model 
structure. The difference between magnetic models for different 
rotational phases is even smaller since g s g varies maximum by a 
factor of S3 2 for Model 1 . 

Finally, we stress that it is difficult to conclude anything 
with certainty regarding the preferable model of the magnetic 
field geometry without additional accurate magnetic observa- 
tions of the 56 Ari. Furthermore, other dynamic processes, such 
as Hall's currents and particle diffusion, may contribute to the 
observed variations of hydrogen lines. These processes can not 
be accounted in our modeling due to their complex nature. 
Nevertheless, similar to the results of the Paper I, in this study we 
demonstrate that the observations can be described with a sim- 
ple geometrical approach under the assumption of strong surface 
electric currents in the atmosphere of a main-sequence mCP star. 

5. Conclusions 

With the use of the high-resolution, phase-resolved observations 
of a magnetic CP star 56 Ari and employing modern model at- 
mosphere technique we detected and investigated variations of 
the Stark-broadened profiles of Ha, HB, and Hy lines in a frame- 
work of a Lorentz force model. Several proofs of the presence of 
significant magnetic pressure in the atmosphere of 56 Ari have 
been found: 

- The characteristic shape of the variation during a full rota- 
tion cy cle of the star corresponds to those described bv lKrollI 
( 1989) and other authors as a result of the impact of a sub- 
stantial Lorentz force (see Paper I and references therein). 

- Numerical calculations of the model atmospheres with in- 
dividual abundances demonstrate that the surface chemical 
spots cannot produce the observed variability in the hydro- 
gen line profiles of the star. 

- Our model shows a reasonable agreement with the obser- 
vations if the outward-directed magnetic force is applied 
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Fig. 6. Magnetic field modulus and Lorentz force parameter as a function of rotation phase for several magnetic field models. 
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Fig. 7. Residual profiles of the Ha, H/3, and Wy lines for each of the observed phase relative to the phase 0.114. Thin solid line - 
observations, theoretical profiles are shown for the following models: thick solid line - B q /B& = -2.5, c%jc\ = 7 model (outward- 
directed Lorentz force), dotted line - B q /Bd = -4.5, c%lc\ = 7 model (outward-directed Lorentz force), dash-dotted line - B q /fid = 
0.5, C2/C1 = 0.1 model (outward-directed Lorentz force), dashed line - B q /Bd = —1.0, C2/C1 = 2 model (inward-directed Lorentz 
force). The residual spectra for consecutive phases are shifted in the vertical direction. The thin dashed line gives the zero level for 
each spectrum. 



assuming the dipole+quadrupole magnetic field configura- 
tion. Unfortunately, due to large uncertainties in the avail- 
able observations of the longitudinal magnetic field it was 
not possible to conclude confidently about the strengths of 
the quadrupolar component. 

Taking into account a variety of possible solutions, we find 
that, to fit the amplitude of a phase-resolved variation in Ha, 
HB, and Hy lines, the magnitude of an induced equatorial 
e.m.f. must be in the range 10~ n - 10~ 10 CGS in case of 



outward-directed Lorentz force and ~ 5 x 10 in case of 
inward-directed one. 



6. Discussion 

56 Ari is the second magnetic CP star for which we detected 
the characteristic variation of the hydrogen Balmer line profiles 
and performed detailed modeling of the Lorentz force effect. Our 
previous target, AOp star 8 Aur, also demonstrated significant 
variability in the hydrogen Balmer lines (see Paper I). However, 
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in the case of Aur, Bor ra & Landstreel (11980b provided a much 
more accurate measurements of the longitudinal magnetic field 
variations, which allowed us to determine more precisely geom- 
etry of its surface magnetic field. Nevertheless, for both stars we 
find that the outward-directed Lorentz force is needed to explain 
observations. However the induced e.m.f, for 56 Ari may dif- 
fer by an order of magnitude, but this remains uncertain due to 
poorly known magnetic field geometry. 

A single-wave variation of the residual spectra is a charac- 
teristic signature for both 9 Aur and 56 Ari. Since both stars 
have high inclination angles and magnetic obliquities, in the 
framework of our Lorentz force model this variation indicates 
the presence of a more complex magnetic field geometry than 
a simple dipole. Such a variation can be obtained in the non- 
dipolar theoretical models by a proper choice of induced e.m.f. 
for each of the multipolar components (e.g. the cilc\ ratio in the 
case of dipole+quadrupole combination). Furthermore, for both 
stars the amplitude of the longitudinal magnetic field variation 
is about w 500 G, which can be the reason for similar amplitude 
of the detected Balmer lines variation (« 1%) since the effec- 
tive temperatures of stars are different (T e ff(9 Aur)= 10400 K, 
r eff (56 Ari)= 12 800 K). 

Similar to the Paper I, we do not consider here any details 
about the physical mechanisms that could be responsible for the 
observed Lorentz force. The final conclusion about the nature 
of the significant magnetic pressure can only be obtained when 
more sophisticated models of the magnetic field evolution and 
its interaction with highly magnetized atmospheric structure will 
become available and/or alternative models will b e tested (how- 
ever, for some of the estimates see discussion in Shulya k et alj 
l2007h . 

In the present work we made use of a simple geometrical 
ID model of Lorentz force: the surface averaged values of the 
transverse magnetic field and the magnetic field modulus are in- 
troduced in the hydrostatic equation of the stellar matter. Future 
investigations can benefit from taking into account 2D effects 
with direct surface integration of the hydrogen line profiles com- 
puted with individual models. This could probably also open a 
possibility to account for the Hall's currents. Unfortunately, as 
it was mentioned above, this is difficult to do at present, but of 
no way impossible once more computational resources become 
available. 

The dependence of the observed variability in hydrogen lines 
upon the magnetic field geometry and strength is one of the key 
element in our investigation. If such a dependence exists, it could 
bring a number of theoretical constrains about the interaction of 
the magnetic field with stellar plasma. So far, we have analyzed 
only two stars with occasionally similar longitudinal magnetic 
field intensity. Note that we are limited to stars for which the con- 
figuration of a magnetic field can be extracted from the literature 
and which can be observed with highly stable spectrometers like 
BOES in order to reduce possible errors in spectra processing. 
Thus, observations of other mCP stars are needed to conclude 
about the connection between magnetic field and variability seen 
in Balmer lines. 
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